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Control-theoretic modeling of the human operator dynamic behavior in manual control 
tasks has a long and rich history. In the last two decades, there has been a renewed interest 
in modeling the human operator. There has also been significant work on techniques used 
to identify the pilot model of a given structure. The purpose of this research is to attempt 
to go beyond pilot identification based on collected experimental data and to develop a 
predictor of pilot behavior. An experiment was conducted to categorize these interactions 
of the pilot with an adaptive controller compensating during control surface failures. A 
general linear in-parameter model structure is used to represent a pilot. Three different 
estimation methods are explored. A gradient descent estimator (GDE), a least squares 
estimator with exponential forgetting (LSEEF), and a least squares estimator with bounded 
gain forgetting (LSEBGF) used the experiment data to predict pilot stick input. Previous 
results have found that the GDE and LSEEF methods are fairly accurate in predicting 
longitudinal stick input from commanded pitch. This paper discusses the accuracy of each 
of the three methods - GDE, LSEEF, and LSEBGF - to predict both pilot longitudinal and 
lateral stick input from the flight director’s commanded pitch and bank attitudes. 


I. Introduction 

S ignificant research has been conducted to model or identify the pilot, as a way to quantify handling 
qualities or to better understand the behavior of a human pilot in controlling a vehicle.^ ^ The classical 
McRuer crossover model of compensatory manual control states that the human operator will adapt to vehicle 
dynamics by providing lead equalization, with an associated computational penalty, such that the combined 
operator- vehicle transfer function is proportional to an integrator at the crossover frequency/ Another 
approach to understanding the human-as-the-controller behavior was the development of the optimal control 
model by Kleinman, Baron, and Levinson based on the assumption that a well-trained and motivated human 
controller behaves optimally in some sense, adjusting the pilot’s compensation for a given vehicle and task, 
subject to human limitations/ 

In the last two decades, there has been a renewed interest in modeling the human operator. Further 
model development by Hess in the frequency domain, a structural model, ^ and Schmidt in the time domain, a 
modified optimal control pilot model, attempted to take into consideration appropriate relevant feedback 
that influenced a pilot’s behavior such as proprioceptive cues, and visual and vestibular feedback. While 
these models proved useful, their application to realistic simulation tasks involved a great deal of complexity. 
In an attempt to simplify the model while still capturing the essential pilot behavioral characteristics, Hess 
proposed a simplified model that has been recently applied to systems with time-varying dynamics. 

In addition to several relevant model structures to describe the pilot-as-a-controller behavior, there has 
also been significant work on techniques used to identify the pilot model of a given structure. These 
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included various applications of parameter identification ranging from wavelets-based approach to frequency- 
based system identification techniques. 

This research is investigating new analytical methods to model the pilots changing behavior over short 
time periods in response to changing aircraft dynamics. This research is attempting to go beyond pilot 
identification based on collected experimental data and to develop a predictor of short-term pilot control 
behavior in real time. The long term goal of this research is to describe a model that goes beyond identifying 
pilot behavior changes due to vehicle dynamic changes but captures changes in pilot’s behavior due to a 
myriad of factors, such as pilot incapacitation. This may enable identification of potential reasons for such 
behavior change. 

This paper is organized as follows: section II describes the experiment, section III gives a brief overview 
of the estimation methods used, section IV discusses the results, and section V offers conclusions. 

II. Adaptive Controller / Pilot Interaction Experiment 

The data this research used is from an experiment that was conducted to categorize the interaction of a 
pilot with an adaptive controller during control surface failures. This experiment is described below. 

The adaptive controller / pilot interaction experiment was conducted to categorize the interactions be- 
tween a pilot and an adaptive controller during a change in aircraft dynamics represented by control surface 
failures. One of the objectives of this experiment was to determine how the adaptation time of the 
controller affects pilots. This was accomplished by considering whether an adaptive controller helps pilots 
during control surface failures (by measuring tracking error) and how this controller affects pilots before, 
during, and after the control surface failures. The control surface failures were either stuck surfaces or slowly 
moving surfaces. A human-in-the-loop experiment looked at the effects of these control surface failures on 
pilot performance during the cruise phase while initiating a climb, descent, or a heading change maneuver. 
These maneuvers were indicated on the primary flight display (PFD) via the flight director and on the en- 
gine instrument displays (EID) horizontal and vertical navigation map displays. Each subject experienced 
four adaptation times: zero seconds, three seconds, seven seconds, and no adaptation (Never). These times 
indicated how long it took the adaptive controller-aircraft system to settle to a new set of dynamics and 
specific adaptation values are based on the response speed of the nominal aircraft dynamics. Zero seconds 
indicated the fastest possible adaptation time, essentially the processor speed. Eor a detailed description of 
this experiment, see Trujillo and Gregory^^’^^ and Trujillo, et al?‘^ 

A. Independent Variables 

The above experiment supplied the data to determine the accuracy of the estimation algorithms in predicting 
subject stick input and the estimation algorithms’ parameter values (estimation gain, forgetting rate, and 
maximum forgetting rate described in Eqs. 7-9 on page 4). 

The primary independent variables were vehicle commanded pitch angle and vehicle commanded bank 
angle. Subjects saw this information displayed on the FED. The difference in ability to predict longitudinal 
and lateral stick was also considered, thus making stick direction an independent variable. 

Subject was also included as an independent variable because it was hypothesized that there might be 
a subject dependency. The above experiment collected flight technical data from 17 subjects, all of which 
were airline transport pilots. 

The above experiment also controlled for time for the adaptive controller to fully adapt. Therefore, 
adaptation time (0 seconds (sec), 3 sec, 7 sec, and Never) was also an independent variable. 

Another independent variable, section, consisted of before, during, and after the failure adaptation. Eor 
an adaptation time of 0 sec, there is no during failure section and for an adaptation time of Never, there is 
no after failure section. 

Table 1 on page 3 details these independent variables. 

B. Dependent Variables 

The estimation methods’ parameter values (PV) are calculated using Eqs. 7-9 below (see Section III). These 
values are the estimation gain from the gradient descent estimator (Eq. 7), the forgetting rate from the least 
squares estimator with exponential forgetting (Eq. 8), and the gain matrix magnitude from the least squares 
estimator with bounded gain forgetting (Eq. 9). The forgetting rate. A, was normalized to a range of 0-1. 
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Table 1. Independent and Dependent Variables 


Independent Variables 

Name 

Values 

Algorithm 

GDE 


LSEEE 


LSEBGE 

Subject 

1-17 

Adaptation Time 

0 sec 


3 sec 


7 sec 


Never 

Section 

Before 


During 


After 

Direction 

Lateral (Gommanded Bank Angle) 


Longitudinal (Gommanded Pitch Angle) 

Dependent Variables 

Name 

Values 

Parameter Value (PV) 

Estimation Gain (P©): 1-100 


Eorgetting Rate (A): 0.01-1.00 


Gain Matrix Magnitude 1-100 


PA during PV before 

PVa6 

PA after PA before 

Average Error 

0-2 

Average Error 

1 Average Error - | Average Error l^g/ore 

Average Error ^6 

1 Average Error |a/ter - | Average Error l^e/ore 


In order to determine the accuracy of the estimation methods, the predicted pilot longitudinal and lateral 
stick input was compared to actual subject longitudinal and lateral stick input collected from the experiment 
described above. Therefore, 

120 

I Predicted Stick Input (t) — Actual Stick Input (t)| 

I Average Error | = (1) 

where t is the data run time (in seconds). | Average Error | was calculated for each run for each subject. The 
stick inputs were calculated for both the lateral and longitudinal directions. Stick movement was normalized 
from —1 to +1. 

Lastly, to see whether there was a discernible difference in the parameter values or the predicted stick 
accuracies during the failure and after the failure as compared to before failure values, these were computed 
in the following manner: 


during ^^bef ore ( 2 ) 

during ^^bef ore ( 3 ) 

I Average Error^f,] = | Average Error^„^i„g - Average Error (4) 

I Average Error^^j = | Average Error^j^^^ — Average Error^g | . (5) 

Table 1 also details these dependent variables. 
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III. Estimation Methods 


This section describes the three estimation methods. The algorithms to estimate pilot longitudinal stick 
input used vehicle commanded pitch attitude data from the experiment described above. The subjects’ 
longitudinal and lateral stick inputs, and longitudinal and lateral flight director commands (commanded 
pitch angle and commanded bank angle respectively) were recorded in the above experiment (section II). 

Consider a pilot input-output model given by a general linear in-parameter form 


y{t) = (6) 

where y(t) is system output (longitudinal or lateral stick in this instance), x{t) is input signal (commanded 
pitch or bank angle in this instance), is the regressor composed of basis functions, and 0(t) is an 

unknown parameter vector. 

By estimating the unknown parameter vector, 0(t), system output is estimated by Eq. 6. Three different 
methods for parameter estimation are considered in building a pilot input-output model. These models are 
a gradient descent estimator (GDE), a least squares estimator with exponential forgetting (LSEEE), and a 
least squares estimator with bounded gain forgetting (LSEBGE).^^ 

The online implementation of the GDE is given by 


kt) 




\ 

\ y{t) ) 


( 7 ) 


where T is the estimation gain. Note longitudinal and lateral stick deflection was normalized to ±1 range. 
Similarly, the online representation of the LSEEE is 


kt) = -r(i)i>(a;(i))(i>(x(i))^0(i)-y(i)), 0(0) = ©o 

t{t) = r(o) = Tq 


The variable A, the forgetting rate, is part of exponential forgetting factor which enables LSEEE to track 
time-varying parameters, albeit slowly varying. The rest of the variables are the same as described above 
for GDE. 

Lastly, the LSEBGE online representation is 

kt) = -T{tMx{t)){^{x{t)Ve{t)-y{t)), 0(0) = ©0 

t{t) = X{t)T{t)-Tmix{t))^T{t), r(0) = To ( 9 ) 

^ q.-M) 

where Aq is the maximum forgetting rate and ko is the pre-specified bound for gain matrix magnitude. The 
rest of the variables are the same as described above for LSEEE. 

Eor GDE, the estimation gain, T, was recursively calculated for each subject’s run until the estimated 
longitudinal and lateral stick inputs matched the actual longitudinal and lateral stick inputs as close as 
reasonably possible. The forgetting rate. A, for LSEEE, was also recursively calculated for each subject’s 
run until the estimated longitudinal and lateral stick inputs matched the actual stick inputs as close as 
reasonably possible. Einally, for the LSEBGE algorithm, the bound for gain matrix magnitude, ko, as 
recursively calculated until a value of A was arrived at that the estimated longitudinal and lateral stick 
inputs reasonably matched the actual stick inputs - once again for each subject. 


IV. Results 

Previous results have found that the GDE and LSEEE methods are fairly accurate in predicting longi- 
tudinal stick input from commanded pitch.^^’^^ These results will discuss the accuracy of each of the three 
methods - GDE, LSEEE, and LSEBGE - to predict both pilot longitudinal and lateral stick input from the 
flight director’s commanded pitch and bank attitudes. The results will also discuss how T, A, and ko vary 
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depending on controller adaptation time during each run’s three segments: before failure, during failure, and 
after failure. 

The above estimators were implemented in MATLAB® and Simulink® R2013b^. Further data analyses 
was done using IBM SPSS® Statistics 22. Significance was set at p < 0.05. 

A. LSEBGF Ao Variance 

As shown above in Eq. 9, A(t) is dependent on Aq, the maximum for- 
getting rate. This value, Aq, was initially set at the average LSEEF 
\{t) calculated from Eq. 8 before the failure occurred, which was 0.49. 

Further analyses indicated that LSEEF A(t) varied by the direction in- 
dependent variable, lateral and longitudinal, before the failure occurred 
(F(i 1963 ) = 12.96; p < 0.01). Table 2 indicates the values Aq was set at when running the LSEBGF algorithm 
again when taking into consideration whether the lateral or longitudinal stick inputs were being computed. 

However, when analyzing PV,PV[^ 5 ^q 5 ], | Average Error | and | Average Error the results for LSEBGF 
with Aq set to 0.49 irrespective of stick direction and for Aq set to 0.52 for the longitudinal direction and 0.46 
for the lateral direction always grouped together when looking at the Tukey Honestly Significance Difference 
(HSD) post-hoc analysis. Therefore, further results using the LSEBGF algorithm will only be reported using 
Ao = 0.49. 

B. Subject Effects 

Not surprisingly, subject was statistically significant for PV and | Average Error |. However, when looking 
at the Tukey HSD post-hoc analysis, there was considerable overlapping in the groups. In general, once 
the overlapping was taken into consideration, subjects appeared to group into two groups - high and low. 
Table 3 details the differences between these two groups. Because the subjects could generally be split up 
into two groups, subject was included in the analysis model for PV and | Average Error | but subject and any 
interactions including subject will be reported using only the two subject groups. 


Table 3. Subject Grouping 




PV 

Average Error 



Group 1 

Group 2 

Group 1 

Group 2 

Section 

Method 

Direction 

Value 

Value 

Value 

Value 

Before 

GDE 


10 

15 

0.0087 

0.0094 


LSEEF 


0.45 

0.53 

0.0311 

0.0515 


LSEBGF 


33 

39 

0.0528 

0.0757 

During 

GDE 


18 

23 

0.0271 

0.0335 


LSEEF 


0.37 

0.47 

0.0766 

0.1016 


LSEBGF 


27 

35 

0.1354 

0.1750 

After 

GDE 


17 

22 

0.0169 

0.0226 


LSEEF 


0.10 

0.10 

0.0664 

0.1124 


LSEBGF 


15 

18 

0.1543 

0.2188 


Table 2. LSEBGF Aq Values 


Section 

Direction 

Ao 

Before 


0.49 

Before 

Longitudinal 

0.52 

Before 

Lateral 

0.46 


When looking at the difference between before and during^ and before and after the failure, subject was 
once again significant but the groupings from the Tukey HSD post-hoc analysis essentially collapsed into a 
single group. Therefore, subject was included in the analysis model for Py[db,ab] | Average Error 
but subject and any interactions including subject will not be reported on. 

C. Ability of Estimation Methods to a Bounded Solution 

The three algorithms were able to predict subject stick input with varying degrees of success. GDE and 
LSEBGF were able to converge to a solution 86% and 76% of the time respectively while LSEEF only was 

^The use of trademarks or names of manufacturers in the report is for accurate reporting and does not constitute an 
official endorsement, either expressed or implied, of such products or manufacturers by the National Aeronautics and Space 
Administration. 
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able to converge to a solution 44% of the time. Convergence was defined as having the predicted output 
within 10% of the actual output over a defined sampling time (section in this case) with a local minimum for 
these differences. The success of the algorithms might have depended on the length of the runs, in this case 
120 sec, and the amount of stick movement during this time, which was not a great deal considering most 
of the 120 sec run was in straight and level flight. Longer sampling times with more stick movement may 
increase the success of the algorithms to converge to a solution. Irrespective of sampling time, the GDE and 
LSEBGF appear to be the best algorithms to use to converge to a solution. 


D. Accuracy of Estimation Methods 

The three estimation algorithms had varying degrees of ability to accurately predict subject stick input. 
This was hinted at in Table 3 on page 5 when looking at the | Average Error | columns. Illustrative examples 
of the algorithms’ ability to predict actual stick commands are shown in Figure 1. 

GDE LSEEF LSEBDF Algortllnni 



(a) GDE Algorithm Longitudinal Stick 
Estimation Example 


(b) LSEEF Algorithm Longitudinal 
Stick Estimation Example 


(c) LSEBGF Algorithm Longitudinal 
Stick Estimation Example 


Figure 1. Representative Longitudinal Stick Estimation Examples with Absolute Error Between Actual and 
Estimated Stick Input 


1. Predicted Stick Input Accuracy Before Failure 

In the short time in the data runs before any failures occurred (<30 sec), subject by algorithm (^( 48 , 10858 ) = 
1.72; p < 0.01) and subject by direction logss) = 2.16; p < 0.01) were significant for | Average Errorj. In 
general, the GDE algorithm was the most accurate in predicting stick input before the failure (Figure 2(a)^ 
on page 7) irrespective of subject. In particular, the GDE method and lateral direction were the least 
sensitive to subject effects (Figures 2(a) and 2(b) respectively on page 7). Thus, if the PV is not specifically 
tuned to the pilot, then the GDE algorithm in the lateral direction provides the most accuracy in predicting 
pilot stick input during normal flight operations; and therefore, the GDE algorithm may be sensitive to 
detecting possible piloting changes during non-normal flight operations. 

2. Predicted Stick Input Accuracy During and After Failure 

As with the predicted stick input accuracy before the failure, the GDE algorithm was generally the most 
accurate, especially during the failure and for lateral stick input prediction (Figures 3(c) and 3(a) respectively 
on 7). Algorithm with adaptation time and direction (^'('g ioi 84 ) = 4.65, p < 0.01), and with adaptation time 
and section (E '('3 ioi 84 ) = 3.46, p < 0.02) (Figure 3 on page 7) were statistically significant. However, the GDE 
algorithm was not the most accurate in predicting longitudinal stick input and after the failure. In this case, 
the LSEBGF method was the most accurate (Figures 3(b) and 3(d) on page 7), albeit not by much. This 
was most likely due to the GDE method having a larger range of absolute average error for some of the data 
runs in the longitudinal direction and after the failure occurred. The errors in the longitudinal direction may 

^The error bars in the graphs indicate one standard error of the mean. 
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GOE LSEEF LSEBGF Lateral Longitudnal 

Algorithm Direction 

(a) Absolute Average Error by Algorithm and Subject (b) Absolute Average Error by Algorithm and Subject 
Group Group 


Figure 2. Predicted Stick Input Accuracy Before Failure 


Lateral Direction 


Longitudinal Direction 



(a) Absolute Average Error by Algorithm and Adaptation 
Time for Lateral Direction 

During Failure 



0 3 7 Never 


Adaptation Time (sec) 

(b) Absolute Average Error by Algorithm and Adaptation 
Time for Longitudinal Direction 

After Failure 



(c) Absolute Average Error by Algorithm and Adaptation 
Time During Failure 



(d) Absolute Average Error by Algorithm and Adaptation 
Time After Failure 


Figure 3. Algorithm Effects on Predicted Stick Input Accuracy 
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be due to the data runs involving an altitude change which persisted for the rest of the data run after the 
altitude change was initiated; therefore, subjects had to deal with the failure and with the altitude change 
for the vast majority of the run. While subjects were constantly adjusting in the longitudinal direction 
for altitude changes, the errors may have been building up in the GDE algorithm whereas the LSEBGF 
algorithm with bounded gain forgetting may have been able to null out earlier large errors. If this is the 
case, then for detecting piloting changes in short runs, the GDE algorithm is still viable; however, for longer 
runs or where there are constant changes in flight path, then LSEBGE may be a better algorithm to estimate 
pilot stick input and to detect changes in piloting because of its forgetting characteristic. 

E. Ability to Detect Changes 

When trying to detect changes in piloting, subject was not significant when analyzing Py[db,ab] | Average 
Error|[^5 Therefore, it may be easier to detect changes in piloting behavior when looking at the differences 
between baseline (in this case, before failure) and later in time (in this case, during and after the failure) 
because individual differences do not play as an important a role. 

1. Absolute Average Error Differenee Between Seetions 

When looking at the absolute average error difference between before, and during and after sections, direction 
by section and adaptation time was significant (E(i 5775) = 7.14,p < 0.01). As seen in Figure 4, it appears 
to be easier to detect changes in the longitudinal direction. When looking at the algorithm by direction 
interaction (E(3 5775) = 4.56, p < 0.01) (Figure 5(a) on page 9) and the algorithm by section interaction 
(^(3,5775) = 20.2, p < 0.01) (Figure 5(b) on page 9), the longitudinal direction and after — before section 
show changes in piloting behavior the most with the LSEBGF algorithm. For the GDE algorithm, the 
direction and section effects are minimal. Therefore, the LSEBGF algorithm may be best for detecting 
changes in piloting behavior when looking at the absolute average error differences between sections in the 
longitudinal direction. 


Lateral Directbn 


Longitudinal Direction 



0 3 7 Never 0 3 7 Never 

Adaptation Time (sec) Adaptation Time (sec) 


(a) Absolute Average Error Difference by Adaptation Time (b) Absolute Average Error Difference by Adaptation Time 
and Section for Lateral Direction and Section for Longitudinal Direction 

Figure 4. Absolute Average Error Difference Interaction 

2. Parameter Value Differenee Between Seetions 

When looking at the differences in parameter values during and after from before section parameter values 
(z.e, during — before and after — before sections), there was an algorithm by direction and section interaction 
(F^s 4615) = 4.89, p < 0.01). As can be seen in Figure 6 on page 9, the GDE algorithm has the most difference 
in parameter values for the during — before and after — before sections. There were also an adaptation 
time and direction, and adaptation time and algorithm interaction with Py[db,ab] (^(2, 5615 ) = 6.63, p < 0.01 
and F(4 5615) = 8.08, p < 0.01 respectively). As seen in Figure 7(a) on page 9, the lateral direction had the 
most change in PV for both during — before and after — before sections. The GDE algorithm also has 
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0.14 





Laltbudn^l Dr^clnn 
1. .-1 LOnUituctinaa Diection 






1 ^ ^ — 

GDE LSEEF LS€BGF 

Algorithm 



LSEEF 

Algorithm 


(a) Absolute Average Error Difference by Algorithm and (b) Absolute Average Error Difference by Algorithm and 
Direction Section 


Figure 5. Algorithm Effects on Absolute Average Error Difference 


Lateral Direction 


Longitudinal Direction 




(a) Parameter Value Difference by Algorithm for Lateral (b) Parameter Value Difference by Algorithm for Longitu- 
Direction dinal Direction 

Figure 6. Algorithm Effects on Parameter Value Difference by Direction 



(a) Parameter Value Difference by Direction and Adapta- 
tion Time 



(b) Parameter Value Difference by Algorithm and Adapta- 
tion Time 


Figure 7. Adaptation Time Effects on Parameter Value Difference 
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fairly large differences in parameter value for all adaptation times (Figure 7(b) on page 9). Therefore, when 
looking at changes in PV, the GDE algorithm in the lateral direction may best indicate changes in piloting 
methods. 


V. Conclusions 

In the last two decades, there has been a renewed interest in modeling the human operator. There has 
also been significant work on techniques used to identify the pilot model of a given structure. The purpose 
of this research is to attempt to go beyond pilot identification based on collected experimental data and to 
develop a predictor of pilot behavior. An experiment was conducted to categorize these interactions on the 
pilot with an adaptive controller during control surface failures. GDE, LSEEE, and LSEBGE algorithms 
used these data to predict pilot stick input. 

The three algorithms were able to predict subject stick input with varying degrees of success. In general 
the GDE and LSEBGE algorithms appear to be the best methods to use to converge to a solution especially 
if the sampling times are short or if there is not much stick movement, such as during straight and level 
flight. 

As for individual differences, the subjects in this experiment appeared to essentially separate into two 
groups - one with lower PVs and | Average Errors | and another with slightly higher PVs and | Average Errors |, 
indicating that some tuning to individuals may be of benefit in order to more accurately predict stick input. 
However, individual differences become a non-factor when comparing the differences between the PV and 
I Average Error | at different times. Therefore, using a standard baseline PV and | Average Error | to compare 
to values acquired during flight may be sufficient to detect piloting changes. Eurthermore, the sampling time 
should be at least 3 and 7 seconds in order to detect a change. If it is shorter than 3 seconds, the change 
may not be sufficient enough to detect. However, if it is too long, the pilot may have time to adapt to the 
conditions such that any differences may be essentially nulled out. 

When looking individually at the prediction of pilot stick input, the GDE algorithm and lateral direction 
movements were the least sensitive to subject effects as long as conditions were not changing significantly and 
the sampling times were short enough such that errors did not build up in the algorithm. Thus, if the PV is 
not specifically tuned to the pilot, then the GDE method in the lateral direction provides the most accuracy 
in predicting pilot stick input during normal flight operations with short sampling times. However, for longer 
runs or where there are constant changes in flight path, then the LSEBGE may be a better algorithm to 
estimate pilot stick input and to detect changes in piloting because of its forgetting characteristic. 

If the goal is to strictly detect piloting changes, then focusing on the difference between baseline values 
and current values should suffice. In this case, the LSEBGE algorithm may be best for detecting changes 
in piloting behavior when looking at the absolute average error differences between sections. However, if 
changes in PV are focused on, then the GDE algorithm in the lateral direction may best indicate changes in 
piloting. 

In summary, the GDE and LSEBGE algorithms have the potential to accurately predict pilot stick 
input from flight director commands. With this capability to predict pilot actions comes the ability to 
detect unexpected changes. Once these changes are detected, then additional data and further analyses may 
determine whether these changes are due to the pilot or a non- normal aircraft. If it is due to a non- normal 
aircraft, then on-board automation may be able to aid the pilot in determining the changing capabilities of 
the aircraft and to try to help alleviate any issues such as by automatically completing electronic checklists. 
If it is due to the pilot, then on-board automation may be able to aid the pilot by changing workload, focusing 
attention appropriately, or, as in the case with an incapacitated pilot, safely land the aircraft. Thus, this 
methodology may lead to a safer national airspace system. 
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